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ABSTRACT 



We present the results of our search for high-redshift Lyman-break galaxies over 
the GOODS-South field. We use HST- ACS data in B, V, i' &l z' , VLT-ISAAC J and 
Qh Ks^ Spitzer-IRAC 3.6, 4.5, 5.8 and S.OiLtm data in conjunction with the new HAWK-I 

Q y-band science verification data to search for dropout galaxies in the redshift range 

^ 6 < z < 9. We survey ^ 119arcmin^ to Yab = 25.7 (5 a), of which 37.5arcmin^ 

^ reaches Yab = 25.9. Candidate z' and Y drop-outs were selected on the basis of a 

colour cut of (y — J)ab > 0-75 mag and {z' — Y)ab > 1.0 mag respectively. We find no 
robust Y-drops {z ^ 9) brighter than Jab < 25.4. In our search for z'-band dropouts 
I {z ^ 6.5 — 7.5), we identify four possible candidates, two with z'-drop colours and 

^ clear Spitzer-IRAC detections and two less likely candidates. We also identify two 

previously-known Galactic T-dwarf stellar contaminants with these colours, and two 
likely transient objects seen in the y-band data. The implications if all or none of 
Psj our candidates are real on the Ultra- Violet galaxy luminosity functions at z > 6.5 are 

explored. We find our number of z'-drop candidates to be insufficient based on the 
expected number of z' drops in a simple no-evolution scenario from the z = 3 Lyman- 
break galaxy luminosity function but we are consistent with the observed luminosity 
function at z 6 (if all our candidates are real). However, if one or both of our 
best z'-dmp candidates are not z > 6.5 galaxies, this would demand evolution of the 
luminosity function at early epochs, in the sense that the number density of ultra- 
violet luminous star- forming galaxies at z > 7 is less than at z ~ 6. We show that 
the future surveys to be conducted with the ESO VISTA telescope over the next five 
^ years will be able to measure the bulk of the luminosity function for both z' and Y 

drop-outs and thus provide the strongest constraints on the level of star-formation 
within the epoch of reionization. 

Key words: galaxies: evolution - galaxies: formation - galaxies: starburst - galaxies: 
high redshift - ultraviolet: galaxies 



1 INTRODUCTION 

The identification of high-redshift galaxies is crucial to de- 
veloping our understanding of the early universe, galaxy 
evolution and the epoch of reionisation. One technique 
for selecting high-redshift galaxies is to search for Lyman 
Break Galaxies (LBGs) through the drop-out technique 
(Guhathakurta, Tyson & Majewski 1990; Steidel, Pettini 
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& Hamilton 1995; Steidel et al. 1996). This relies on find- 
ing a significant flux decrement between two broadband 
fllters indicating a spectral break short-ward of Lyman-a 
at I215.7A in the rest frame of the galaxy, attributable to 
Lyman-a absorption from intervening clouds of neutral hy- 
drogen at lower redshifts (the Lyman-a forest). In higher- 
redshift galaxies this "Lyman break" is pushed to longer 
wavelengths. 

Recent years have witnessed two key developments in 
our understanding of the initial phases of galaxy formation. 
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The first is tiie extension of tfie Lyman-break selection tecfi- 
nique to fiigfier redsfiifts, so tfiat broad-band colours can 
be used to identify galaxies at z = 6 (Stanway, Bunker & 
McMahon 2003; Bunker et al. 2004). The Lyman-break 
drop" selection has been spectroscopically confirmed to find 
galaxies at z ^ 6 (Bunker et al. 2003; Stanway et al. 2007). 
The second key development is the revelation that z > 6 is a 
crucial transition epoch in the history of the Universe. The 
discovery of complete absorption by the Lyman-a forest (the 
Gunn-Peterson trough) in the spectra of quasi-stellar objects 
(QSOs) at z > 6.2 (Becker et al. 2001, Fan et al. 2002) indi- 
cates a large neutral fraction of hydrogen, and the Wilkinson 
Microwave Anisotropy Probe (WMAP) results suggest that 
the Universe was entirely neutral at 2; > 10 (Kogut et al. 
2003) with a protracted period of reionization (Dunkley et 
al. 2009). Hence, exploring the epoch 6 < ^ < 10 is crucial 
if we are to understand what reionized the Universe, and 
thus set the stage for galaxy formation at the end of the 
"Dark Ages". Specifically, ultra-violet (UV) light from star 
forming galaxies during this era has been proposed as the 
most likely reionization mechanism, as the number density 
of high-redshift active galactic nuclei ( AGN) appears too low 
to be solely responsible (e.g. Dijkstra et al. 2004). 

Only a handful of galaxies are known at z > 6.5 
(Bouwens et al. 2005, 2008; Ota et al. 2008), selected through 
Lyman-break broad-band imaging and also narrow-band 
imaging for Lyman-a line emitting galaxies. The highest 
redshift spectroscopically confirmed galaxy currently known 
was discovered by lye et al. (2006) at z = 6.96 and the high- 
est redshift object is the gamma-ray burst found at z ~ 8.3 
(Tanvir et al. 2009). The frontier in spectroscopically con- 
firmed galaxies and large robust samples of Lyman Break 
galaxies is currently z ~ 6. The data from the new VLT 
HAWK-I camera (Pirard et al., 2004; Casah et al., 2006) 
in the y-band filter in conjuction with the availability of 
ISAAC images over most of the ACS/WFC GOODS-South 
area enables us to push this work even further in lookback 
time and to explore the redshift range 6.5 < z < 9. Con- 
straints have also been placed on the Ultra Violet luminosity 
function at ^ 7 by Mannucci et al. (2007). 

To explore this population further we need a large sam- 
ple of these galaxies. This has been hampered until now by 
a lack of sensitivity in the infra-red and the small fields of 
view available. Now with HAWK-I, an instrument with a 
large field of view (7.5 x 7.5 arcmin^) on the VLT, an 8- 
metre class telescope and critically with the y-band filter, 
we can begin to increase the number of Lyman Break candi- 
dates at z > 6.5. The F-band filter is particularly important 
because it is a good discriminant of spectral breaks (a sharp 
cutoff) as opposed to dust reddening, which appears as a 
more gradual slope, as it is centered on l/xm which is close 
in wavelength to the z' filter's peak transmisison wavelength 
of ~ 0.9/im (see Figure[T]). The y-band filter is also currently 
being used by UKIRT Infrared Deep Sky Survey (UKIDSS; 
Lawrence et al. 2007) to search for QSOs at similar epochs 
(e.g. Venemans et al. 2007). 

The structure of the paper is as follows. In Section [2] 
we describe the imaging data. The construction of our cata- 
logues, survey completeness, our z' and y-drop selection and 
the likely contaminants are discussed in Section [3] In Sec- 
tion |4] we discuss our candidates and infer constraints on the 
rest-frame UV luminosity function of star- forming sources. 



Our conclusions are presented in Section [5] Throughout we 
adopt the standard "concordance" cosmology of Qm = 0.3, 
Qa = 0.7, and use /170 = ifo/70kms"^ Mpc"^ We use the 
AB magnitude system throughout (Oke & Gunn 1983). 



2 OBSERVATIONS AND DATA REDUCTION 

We analysed y-band observations from the ESO/VLT 
archive obtained as part of HAWK-I science verification 
program 60.A-9284(B) (Fontana et al. and Venemans et 
al. - "A deep infrared view on galaxies in the early Uni- 
verse"). Two areas within the GOODS-South field were im- 
aged, with centres 03:32:41.0 -27:51:45 (pointing 1) and 
03:32:29.6 -27:44:37 (pointing 2, both J2000). The seeing 
in y-band was in the range 0'.'4 - 0V7 (FWHM) in the in- 
dividual exposures, with the stacked image having seeing of 
^ 0%. 

Individual exposure times were 30 seconds, with 10 such 
exposures averaged to form a single co-added frame of 
300 seconds. A typical sequence was 12 exposures (Ihour), 
although we used sequences with between 6 and 13 expo- 
sures. The telescope was dithered by 5-10^^ between expo- 
sures. There were 18 sequences in pointing 1 (comprising 195 
frames in all, a total of 16.25 hours), and 14 towards pointing 
2 (comprising 138 frames, amounting to 11.5 hours). 

The HAWK-I camera was read out in 4 quadrants, each 
of 2048 X 2048 pixels with a scale 0^.^106 pixel" \ and initially 
these quadrants were reduced individually before final mo- 
saiking. We subtracted the dark current and flat-fielded the 
data using the average of y-band twilight sky flats from 
the ESO HAWK-I archive. For each sequence of ^^12 ex- 
posures, we then used the DIMSUM package within IRAF to 
background-subtract the sky, using the xmosaic task. We 
measured the image shifts interactively using a number of 
compact, bright but unsaturated sources. Where there were 
12 frames in a sequence, we used the average of the 5 frames 
before and 5 after (with objects masked) to determine the 
sky level. DIMSUM also rejected many of the cosmic rays 
present in the images. However, there were various electronic 
read-out artifacts visible which were pernicious and not re- 
moved by DIMSUM. Hence when we combined the 14-18 
final frames output by DIMSUM, we used IRAF.imcombine, 
weighting by the exposure time, with the ccdclip rejection 
algorithm, using the gain and readnoise properties of the 
detector and the Poisson noise of the sky background. A de- 
tector artifact which was not effectively eliminated by this 
rejection was a cross-talk effect, whereby ghost images ap- 
peared in the same detector row as bright objects. As the 
apparent spatial position of the cross talk artifacts on the 
sky remains fixed with respect to the objects (i.e. it moves 
with the objects on the detector during the dithering pat- 
tern), this was not rejected in the co-addition. We eliminated 
this later by visual inspection and by masking detector rows 
affected by bright objects. 

For each quadrant and for both pointings we then de- 
termined the astrometry through comparison with the ESO 
ISAAC J-band image (Retzlaff et al. in prep. )^usmg 60-110 

^ We used version 2 - available from 
http : / / archive, eso . org / archive / adp /GO CD S /ISAAC .imaging _v2 .0 
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objects in each quadrant. A quadratic fit to the distortion 
produced residuals of ^ 0'!l between the coordinates in the 
J- and y-band images. The iraf . wregister routine was 
used to hnearly interpolate the y-band images to map on 
to the same pixel grid as the J-band, with a pixel scale of 
0'.'15. 

We summed the exposure maps output by DIMSUM 
using the same measured shifts, creating a map of the to- 
tal exposure time as a function of position on the sky, and 
correcting the astrometric distortions as described above. 
This was then used to inverse-variance weight the images 
when we combined the four quadrants in both pointings 
to form a final image mosaic. The 119arcmin^ of the Y- 
band image covered most of the ESO J-band and HST ACS 
GOODS images (97.5% of the HAWK-I F-band image over- 
lapped with the ESO J-band). We used the GOODS team 
reductions of the ACS images (Giavalisco et al. 2004), con- 
sisting of F450W B-band, F606W V-band, F775W i'-band 
and F850LP z-hand. The GOODS images had been drizzled 
from the original ACS pixel scale of 0V05 on to a grid of 0^.^03 
pixels. We used version 2.0 of the ACS GOODS image^and 
adopted the AB magnitude zeropoints re-determined for the 
v2 release. 

At the time of writing, the formal y-band zero point for 
HAWK-I was not available, so we determined the photomet- 
ric zero point of the y-band by measuring the Y — J colours 
of objects in identical apertures of 1 arcsec diameter, and 
setting the average AB colour of sources with flat spectra 
between the z' and J-band to be zero. Using our computed 
zero point, the AB magnitude in the y-band is given by 

Yab = 30.46 — (ap. corr.) — 2.5 log^Lol^^unt rate) mag 

where "ap corr" is the aperture correction in magnitudes, 
and "count rate" is the number of counts per second 
recorded. We determined the aperture correction using a 
1 arcsec-diameter aperture to be 0.4 mag for compact but 
unsaturated sources (see Section 3.1). 

The 5cr detection limit for a compact source in a 
1 arcsec-diameter aperture is Yab = 25.7 mag, however 
pointing 1 is slightly deeper and reaches Yab — 25.9 mag 
(5cr). In all, we survey 37.5 arcmin^ to a maximum 5 a depth 
of Yab — 25.9, and an area of 90.6arcmin^ to Yab < 25.7. 
The total area surveyed to Yab < 25.5 was 115.6 arcmin^. 

In order to identify y-band and ^^-band drop-outs 
through their extreme (Y — J) or {z' — Y) colours, we com- 
pare the HAWK-I y-band image with the GOODS team 
J-band and 2;^-band images (taken with VLT-ISAAC and 
HST- ACS respectively). The 5 a limiting magnitudes in 
GOODS ACSv2 measured in l^.^O diameter apertures are 
Bab = 27.2, Vab = 26.90, Zab = 26.09 and Zab = 
26.14. We block-averaged the drizzled 0'.'03 ACS pixels 
5 X 5 to produce a z'-hand frame which was registered to 
the VLT/ISAAC J-band. The limiting magnitudes in I'.'O 
apertures for the VLT-ISAAC images are J = 25.2 and 
Ks = 24.7 (5(7 AB magnitudes). We also made use of 
imaging of the GOODS-South field obtained with the In- 
frared Array Camera (IRAC; Fazio et al. 2004) on board the 
Spitzer Space Telescope which was conducted as part of the 



2 The GOODS ACS v2 images are available at 
|http : //archiv e . stsci . edu/pub/hlsp/goods/v2/j 
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Figure 1. The ACS z' - and HAWK-I y-band passes over plotted 
on the spectrum of a 2; = 6.8 galaxy with 100 Myr of constant 
star formation and 0.2 Solar metallicity (solid line), illustrating 
the utility of our two-filter technique for locating z ^ 7 sources. 

GOODS Legacy program (PID 194, PI Dickinson). IRAC 
uses four broad-band filters with central wavelengths at ap- 
proximately 3.6/im, 4.5/im, 5.8/im and 8.0/im (channels 1- 
4). The data were taken over two observing epochs, with the 
telescope roll angle differing by 180°, and we used the v2 and 
v3 reductions of the GOODS team, with the data drizzled 
onto a 0^.^6 grid from the original 1^2 pixels. We measure 5cr 
limiting magnitudes of 24.76, and 24.87 for IRAC channels 1 
and 2 respectively, measured in 2''.^ 4 diameter apertures, and 
22.77 and 22.81 for IRAC channels 3 & 4 measured in 3'.'0 
& 3^.^ 7 diameter apertures respectively (these limits include 
aperture corrections of ^ 0.7 mag appropriate for unresolved 
sources, e.g. Eyles et al. 2005; 2007). 



3 SELECTION OF Z -DROP AND y-DROP 
CANDIDATES 

3.1 Construction of Catalogues 

Candidate selection for all objects in the field was performed 
using version 2.4.6 of the SExtractor photometry package 
(Bertin & Arnouts 1996). As we are searching specifically 
for objects which drop-out at short wavelengths through 
the Lyman-a forest absorption, we used SExtractor in dual- 
image mode, detecting objects in the longer- wavelength 
band and measuring the photometry within the same spatial 
apertures in the drop-out band(s). We produced separate 
catalogues for the F-band drop outs (using the J-band as 
the detection image) and the z'-hand drop-outs (using the 
y-band as the detection image). To reduce the number of 
spurious sources in the noisy edge regions (where few frames 
overlap) we used the exposure maps as input weight maps 
for SExtractor. For object identification, we adopted a limit 
of at least 5 contiguous pixels above a threshold of 2a per 
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pixel (on the data drizzled to a scale of 0Vl5 pixel"^). Spuri- 
ous detections close to the noise limit were later eliminated 
through colour cuts and visual inspection. 

As high redshift galaxies in the rest-UV are known to be 
compact (e.g., Ferguson et al. 2004; Bouwens et al. 2004), we 
used fixed circular apertures l^.^O in diameter to select our 
candidates and corrected the aperture magnitudes to ap- 
proximate total magnitudes for each filter through an aper- 
ture correction, determined from bright compact sources. 
These were measured to be 0.07 mag in z^-band, 0.4 mag in 
y-band and 0.4 mag in J-band. For compact sources, this 
approach is simpler and more reproducible than using the 
SExtractor curve-of-growth total magnitudes. 

Section 9 in the version 2.0 of the ISAAC J-band image 
produced a number of spurious sources when running SEx- 
tractor in dual image mode and weighting on the exposure 
map. This appeared to be due to an overconfidence in the 
weight map on section 9. This was corrected for by weight- 
ing on the background for this section alone and using the 
weight map for the rest of the J-band image. 
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3.2 Completeness 

The completeness corrections for both the Y— and the J- 
band images were measured in the following way. Approx- 
imately 5000 artificial compact sources were created with 
diameters of 3 pixels and spanning magnitudes between 20 
and 30. These objects were then convolved with the PSF 
and added in to the original image. The new images were 
run through SExtractor again, using the same criteria that 
was employed to generate our object list. The resulting cat- 
alogues were compared with the list of input 'fake ' sources 
and a detection was considered to be made if a source was 
found within 5 pixels of its input position and had a mag- 
nitude correct to within a factor of 2 of the input flux. The 
number of detected sources was then compared to the num- 
ber of input sources for each magnitude bin. 

As described earlier, the F-band image consisted of two 
individual pointings of unequal depth, to measure an accu- 
rate completeness limit, the calculations were determined 
for both pointings separately. The filled circles in Figure [2] 
show the y-band image for Pointing 1 is 95 per cent com- 
plete down to a magnitude of Yab = 24.0 and is 50 per cent 
complete at a magnitude of Yab = 25.9 over the deepest 
area. The l^-band image for Pointing 2 (denoted by open 
diamonds in figure [2| is ^ 95 per cent complete down to a 
magnitude of Yab = 24.0 and is 50 per cent complete at 
a magnitude of Yab = 25.7 over the deepest area in that 
pointing. We also estimate the completeness of the J— band 
image with a similar method. Figure [3] shows the J-band 
image is ~ 90 per cent complete to Jab = 24 and is 50 per 
cent complete at Jab = 25.4. 

We calculated our 5 a limit in the F-band to be 25.9 in 
pointing 1, which agrees with our 50 per cent completeness 
limit (Fig. |2|. We tested the SExtractor parameters by re- 
ducing our detection thresholds and found they had little or 
no effect on our completeness implying that our original pa- 
rameters of at least 5 pixels at 2 a or above were reasonable 
and did not eliminate credible sources from our selection. 



Figure 2. The y-band completeness with SExtractor parameters 
of at least 5 pixels with S/N > 2cr. Pointing 1 is denoted by the 
filled circles and Pointing 2 by the open diamonds. 
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Figure 3. The J-band completeness with SExtractor parameters 
of at least 5 pixels with S/N > 2a 

3.3 z-drop Candidate Selection 

For our 2;^-drop selection, our SExtractor catalogues revealed 
278 candidates with aperture corrected colours {z' — Y)ab > 
1 mag and S/N > 5 in the Y band. We expect that many 
of these candidates will be low-redshift interlopers such as 
low-mass stars or red galaxies at z ^ 2, which can pro- 
duce large (z' — Y) colours (in particular due to 4000 A and 
Balmer breaks - see Figured. We chose this colour cut of 
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Table 1. Properties of our 4 z^-drop candidates, magnitudes are listed with an aperture correction applied as described in the text. 
Where the candidate is undetected we have quoted the 2cr limiting magnitude. Object 2200 is probably a lower redshift contaminant at 
3.25 < z < 3.85 as discussed in the text. Objects 9136 and 9697 are our most convincing 2;^-drop candidates. 



Our ID 


2200 


9136 


9266 


9697 


RA & Dec 


03 32 25.3 -27 52 30.7 


03 32 17.4 -27 43 43.0 


03 32 19.2 -27 43 33.4 


03 32 22.7 -27 43 00.8 




26.87 ± 0.31 


27.19 ± 0.38 


>28.1 


27.84 ± 0.65 


Yab 


25.70 ± 0.14 


25.90 ± 0.18 


25.94 ± 0.19 


25.29 ± 0.10 


Jab 


25.23 ± 0.16 


24.98 ± 0.23 


>26.18 


26.0 ± 0.42 


Kab 


23.8 ± 0.07 


24.98 ± 0.32 


> 25.42 


>25.42 


3.6ab 


21.4 ± 0.01 


23.6 ± 0.054 


>25.76 


24.01 ± 0.07 


^•^AB 


21.1 ± 0.01 


23.55 ± 0.07 


>25.87 


24.35 ± 0.13 


5.SAB 


20.53 ± 0.02 


>23.77 


>23.77 


> 23.77 


S.Oab 


20.46 ± 0.02 


22.24 ± 0.28 


>23.81 


>23.81 


{z' - Y)ab 


1.17 


1.29 


2.10 


>2.47 




redshift 



Figure 4. Model colour- redshift tracks for galaxies with non- 
evolving stellar populations (from Coleman, Wu &; Weedman 1980 
template spectra). The contaminating 'hump' in the {z' — Y) 
colour at z ~ 2 arises when the Balmer break and/or the 4000 A 
break redshifts beyond the i^-filter. 

{z' — Y)ab > Imag because although it does not elimi- 
nate all of the low redshift interlopers it does omit a sig- 
nificant fraction without excluding potential candidates at 
z > 6.5. We show the distribution z' — Y colours for all of 
our detected objects in Figure |5] In order to eliminate ob- 
vious low-redshift contaminants, we compared our list of z- 
drop candidates to the GOODS MUSIC catalogue (Grazian 
et al. 2006) with a matching radius of 0^.^36. This cata- 
logue provides photometry from HST-ACS, Spitzer-IRAC , 
and ground-based [/-band and VLT-ISAAC JHKg imag- 
ing, with PSF-matching used to determine accurate colours. 
The GOODS-MUSIC catalogue includes photometric red- 
shift estimates derived from the 14-band photometry, and 
the catalogue is a combination of a iCs-band and 2;^-band 
selection. 

Of the 278 candidates, 101 appeared in the GOODS- 
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Figure 5. z' — Y colours of all objects in the HAWK-I data. The 
plus signs are 2 a lower limits on the {z' — Y) colour where objects 
are undetected in the z^-band. 

MUSIC catalogue, mostly with photometric redshifts 
Zphot = 1 — 2.5, although there were two with ZpUot — 6.9 
which were identified by Mannucci et al. (2007) as brown 
dwarfs - these objects are further discussed in Section [43] 

There was also one other candidate with Zphot > 5 that 
was identified in our original search but was found to have 
a match to GOODS-MUSIC catalogue object 30046. It has 
detections in the i\ z\ Y, J, Ks and Spitzer hands with very 
strong emission at the Spitzer wavelengths. Its photometric 
redshift of Zphot = 5.14, along with its detection in the i- 
band and its non-detection at the shorter ACS wavelengths, 
leads us to classify this object as a possible F-drop and 
therefore remove it from our sample of z'-drop candidates. 

We confirmed that all of the Zphot < 5 matches to 
GOODS-MUSIC had detections in one or more of the deep 
ACS B-, V- and i^-bands, ruling out high-redshift interpre- 
tations due to the absence of a break at Lyman-a. In Fig- 
ure M we show the {z' — J) colours of our z'-drop candi- 
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Figure 6. z' — J colour versus the photometric redshift derived 
from the GOODS-MUSIC catalogue. The z^-drop out candidates 
that satisfy our colour selection criteria are marked with crosses. 
All of the sources in the GOODS-MUSIC catalogue with S/N > 3 
in the J-band are denoted by dots 



dates with GOODS-MUSIC matches over-plotted with the 
full GOODS-MUSIC catalogue (with a J-band threshold 
of S/N > 3). As can be seen, most of our sources have 
{z' — J) MUSIC > 1-0, as would be expected from our selec- 
tion of {z' — Y) > 1.0, and the bulk lie at Zphot ~ 2 (as 
would be expected for the interloper populations - see Fig- 
ure [i]). The 2;^— drop candidates that do not lie at 2; ^ 2 
or z > 6 are largely attributable to the GOODS-MUSIC 
catalogue dealing with total magnitudes, whereas we chose 
to use aperture magnitudes (more accurate for the expected 
compact nature of high-redshift galaxies), with large low- 
redshift galaxies having a greater aperture correction than 
we adopted. Also, colour gradients within galaxies mean that 
aperture photometry may select red regions of galaxies (e.g. 
spiral bulges) as z'-drop candidates; the HST-ACS z^-band 
has better resolution than the ground-based F-band so we 
are sometimes susceptible to edges of large objects as spuri- 
ous candidates. 

For the 177 z^-drop candidates which did not lie within 
0'!36 of a GOODS-MUSIC source, we visually inspected all 
four HST-ACS wavebands, as well as the HAWK-I F-band 
and the ISAAC J and Kg bands, to ascertain whether the 
y-band detection was real, and if there was any detection at 
other wavelengths. Flux in the ACS B-, V- or i^-bands would 
be incompatible with the source being a z^-drop Lyman- 
break galaxy at z > 6.5. We found that 17.5 per cent of 
the remaining 177 candidates were detector artifacts (most 
frequently the cross-talk effect due to a bright object in the 
same detector row, manifesting a positive-negative dipole 
signal). Ghost image haloes around bright stars accounted 
for another 14 per cent, and 14 per cent again of the can- 
didates were unreliable due to falling in regions of excess 



noise (despite using exposure weight maps to cut-down on 
spurious detections). 

To verify the reality of our six remaining z-diop can- 
didates, we split the Y band data in to two halves (in time) 
and combined the first half and the second half of the data 
separately. We then measured the Y band magnitudes of 
the candidates in both halves of the data (along with some 
reference star (in case the seeing or magnitude zeropoints 
differed over time) to check for consistency. Two of the six 
remaining candidates were eliminated during this process as 
they were only visible in the second half of the data. These 
objects appeared very bright in the F-band and were un- 
detected in all of the other bands. This prompted us to ex- 
amine the individual images which revealed the two sources 
to be time variable. Each night's data was combined sep- 
arately and the photometry on the individual nights were 
measured for each candidate. Both objects were undetected 
on 2007-10-17, 2007-10-18, and 2007-10-19 but were visi- 
ble on 2007-12-01 and 2007-12-02 with magnitudes from 
Yab = 24.77 — 24.88. The two transient objects have co- 
ordinates of a = 03^32^54.4" 6 = -27^53^35.7" and 
a 03^32^33.5" 6 -27^49^38.3". 

Hence we find that 46.7 per cent of the 177 z'-drop 
candidates without GOODS-MUSIC matches are spurious. 
A comparable fraction (50 per cent) had detections visi- 
ble in the HST-ACS images; it is probable that they did 
not have corresponding GOODS-MUSIC matches because 
of the small matching radius we adopted (0^.^36) to cut down 
on multiple matches. Small astrometric shifts in some re- 
gions, the 2;^-band and X-band magnitude cuts in GOODS- 
MUSIC, and spatially extended galaxies with colour gradi- 
ents account for these objects not having GOODS-MUSIC 
matches. Hence, from an initial colour selection resulting in 
278 objects, we have 4 remaining z^— drop-out candidates 
without GOODS-MUSIC matches. Table [l] contains our list 
of z^-drop candidates. 

Some of our z^-drop candidates may simply meet our 
selection criteria due to photometric scatter. To assess how 
significant this may be in our dataset we use our parent cat- 
alogue of y— band detected sources over the 90.6 arcmin^, 
which reaches a 5cr depth of Yab = 25.7 mag (see section |2|, 
and randomly redistributed the magnitudes according their 
corresponding uncertainties in both the z' and y— bands 
and calculated the new z' — Y colour for each source. We 
find, on average, that we detect nine objects purely due to 
photometric scatter in the parent catalogue. However, if we 
then take the fraction of the total sources that we exclude 
through our subsequent cuts through matching to the HST- 
ACS images etc., then we expect approximately 0.2 sources 
to be spurious in our final candidate list. Therefore, it is con- 
ceivable that one of our sources is a product of photometric 
scatter. Due to the lack of detection in any band other than 
the y-band then Object 9266 is plausibly spurious. 

3.4 y-drop Candidate Selection 

The y-drop candidate selection was carried out in the same 
manner as the z^-drop selection. Our criteria for the y-drop 
candidates were a colour difference of (Y — J)ab > 0.75, 
a signal-to- noise ratio 5'/A^ > 5 in the J-band and a value 
in both the Y- and J-band exposure maps equivalent to a 
minimum of 2.5 hours of observation in the Y band. This se- 
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Table 2. Properties of our 4 y-drop candidates, magnitudes are listed with an aperture correction applied as described in the text. 
Where the candidates are undetected in Y we have quoted our 2a limiting magnitude 



Our ID 


2058 


4551 


5512 


4532 


RA & Dec 


03 32 27.6 -27 51 04.1 


03 32 16.2 -27 47 39.1 


03 32 27.5 -27 46 14.7 


03 32 48.3 -27 47 39.9 


Yab 


>26.6 


>26.6 


>26.6 


27.47 ± 0.6438 


Jab 


25.37 ± 0.19 


25.07 ± 0.19 


25.05 ± 0.19 


25.37 ± 0.19 


Kab 


>25.42 


>25.42 


>25.42 


>25.42 


3.6as 


>25.76 


>25.76 


>25.76 


>25.76 




>25.87 


>25.87 


>25.87 


>25.87 


5.8as 


>23.77 


>23.77 


>23.77 


> 23.77 




>23.81 


>23.81 


>23.81 


>23.81 


{Y - J)ab 


>1.23 


>1.53 


>1.55 


2.1 



lection yielded a list of 133 possible y-drop candidates. We 
then compared our list to the GOODS MUSIC catalogue to 
eliminate the low redshift interlopers from our selection. In 
all, 98 of our 133 objects had matches within 0^.^36 of our 
candidates. One of these objects had no detection in the 
bands V oi i' so we kept it in our candidate list. This 
resulted in 37 y-drop candidates. These remaining objects 
were inspected more closely with postage stamps in B, 
i' ^ z' ^ y, J and Kg bands. We found 16 per cent of our re- 
maining candidates to be ghost image haloes around bright 
stars, 27 per cent were detections picking up on the edges 
of bright galaxies in the J-band. Another 3 per cent fell on 
noisy regions of the y-band and 43 per cent of our candi- 
dates had visible ACS detections. From our original list of 
133 candidates only 4 possible y-drops remain, and these 
are listed in Tabled 



4 DISCUSSION 
4.1 Candidates 

We have four z^-drop candidates remaining after eliminating 
artifacts and low-redshift interlopers. Our candidates span 
a range Yab = 25.3 — 26.0 (after applying our aperture cor- 
rection), two of which have > 2cr detections in the ^^-band. 
Three of the candidates have strong detections in the IRAC 
wavebands. We obtained photometric redshift s for these ob- 
jects because the detection at the IRAC wavelengths in- 
creases the likelihood of a valid detection and also improves 
the accuracy of the photometric redshift solutions. We used 
the publicly available software Hyperz (Bolzonella, Mi- 
ralles & Pello et al. 2000) to derive our photometric redshift 
estimates in the redshift range z = — 9. We allowed vi- 
sual extinction values between Ay =0 — 4 and assumed the 
Calzetti (1997) reddening law. We used the eight Bruzual & 
Chariot (2003) template spectra provided, using solar metal- 
licity and included the following 11 filters V, i\ z\ y, 
J, Ks, and the four IRAC channels . We employed option 
2 for the error treatment of an undetected source which as- 
sumes that the flux in that filter and its la error are equal 
to half the flux of the limiting magnitude (i.e. the error bar 
ranges from flux=0 to the la limiting magnitude in that 
waveband). 
Object 2200 



Hyperz is available at http://webast.ast.obs-mip.fr/hyperz/ 



This object displays a strong detection in the Y—, J— 
and Ks— hands as well as a significant detection in the IRAC 
channels 1 and 2. There is some flux detected in the z'- 
band, however this is to be expected for some candidates 
as an examination of Figure [l] shows the z— and y— filter 
transmission curves do overlap significantly. This means as 
the Lyman break moves through the y— band filter with in- 
creasing redshift its contribution to the 2;^-band flux will de- 
crease but may not entirely disappear. The Ks— hand source 
is slightly offset from the detection in the other wavebands 
for object 2200, this prompted us to widen our search area in 
the GOODS-MUSIC catalogue to a I'.'O radius. This larger 
radius yielded a match to an object in GOODS-MUSIC with 
an ID 30199 and Zphot = 2.73. This object was also identified 
by Stanway et al. (2008) as a possible z^-drop. However, the 
detection by Stanway et al. is centered 0^.^8 from our source. 

This object also has a reported MIPS 24 /xm detection 
(source mip003485 in Alonso-Herrero et al. 2006), and coin- 
cides with a 0.5 — 2keV Chandra X-ray source, and hence 
probably has an AGN contribution. This object has also 
been presented in Dunlop, Cirasuolo & McLure (2006), their 
object 2336, who derive a photometric redshift in the range 
z = 3.25 — 3.85. The photometry presented in Table [l] is for 
an object at the position of the y — band source, this is offset 
from the i^s- band source and as such has a fainter magni- 
tude than that given by Dunlop et al. (2006). We believe 
that the y— band detection is probably associated with the 
iT- source, although this requires additional imaging and/or 
spectroscopy to confirm. 
Object 9136 

This candidate displays a strong detection in the Y 
band and is also detected in the J and i^s bands, and like ob- 
ject 2200 it is strongly detected at the IRAC wavelengths. 
In Figure [7| we show the photometric data points for this 
object with the best fit galaxy template overlaid. The best 
fit solution from Hyperz is Zphot = 7.01 with a secondary 
peak in the probability distribution at Zphot = 7.23 (Fig.|8|. 
This best fit solution is for a starburst galaxy with Ay = 0. 
We calculated star formation rates for our objects based on 
the rest frame UV continuum at 1500 A (see Bunker et al. 
2004). For object 9136 we find a star formation rate of ~ 22 
Moyr-^ 
Object 9266 

This object is detected in the y-band, but not in the J- 
and i^s -bands. However the limits in these bands are fainter 
than the measured y-band magnitude. This could indicate 
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Figure 7. Best fit SED at Zphot = 7.01 to Object 9136 with 
photometry overlaid and 2 sigma upper limits denoted by down 
arrows 




redshift 

Figure 8. The redshift probability distribution for Object 9136, 
showing the best-fit photometric redshift of z = 7.01. 

that the object is spurious or a result of line contamina- 
tion in the y— band filter, or simply that its continuum is 
fainter than the J— and Ks—hs^nd limits but bright enough 
to be detected in the F-band. We did not fit a photometric 
redshift for this object due to its limited detections. 
Object 9697 

This object was previously identified as an i^-band drop- 
out in the GOODS ACSvl data by Bouwens et al. (2006), 
and is #2226643007 in their catalogue, with z'j^b = 27.54 ± 



23 
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Figure 9. The best-fit SED at Zphot = 6.92 for Object 9697 
with the multiband photometry overlaid. The 2cr upper limits are 
denoted by the arrows 

0.18, [i — z') > 1.3 and infrared magnitudes Jab — 26.04 
and Kab > 25.4. The source is detected in the Y- and J- 
bands and there is also a strong Spitzer detection. However 
there is a nearby source in the z^-band unconfirmed at the 
other wavelengths and another detection ^ 0V7 to the east 
in the B and V bands which, while unassociated with the z- 
band detection, may be at least partially responsible for the 
IRAC flux as it does fall within the IRAC aperture. However 
it would be a rather unusual object to be detected in B and 
V and the Spitzer hands and undetected in z^, y, J, and Kg. 
The i' z'Y JKs colours appear consistent with a high-redshift 
interpretation. We also ran Hyperz on this object assum- 
ing the Spitzer flux was contributed by our candidate to 
determine its photometric redshift. In Figure [9] we show the 
photometry of this object with the best-fit galaxy template 
overlaid. The best-fit solution was ZpUot — 6.92 with a sec- 
ondary peak in the probability distribution at Zphot = 5.22, 
see Figure [To] This best-fit solution is for an elliptical galaxy 
with Av — 0. We calculate the star formation rate based on 
the UV continuum to be ^ 25M0 yr~^. 

Another possibility, other than a high redshift interpre- 
tation, is that the objects which are strongly detected in the 
IRAC bands, 9136, 9697 and 2200 may be similar to IRAC- 
selected extremely red objects (lEROs) see (Yan et al. 2004). 
But in the sample discussed by Yan et al. the sources had 
optical detections which our candidates do not. At the high- 
redshifts that we estimate for these galaxies, if there was an 
old stellar population present, the 4000 A break could fall 
between the Kg and the 3.6/im bands. This could explain 
some of our sources increasing in brightness at the IRAC 
wavelengths. 

We also explored the possibility that our remaining can- 
didates could be brown dwarfs. Patten et al. (2006) pre- 
sented observationally-derived colours for various spectral 
types of M, L and T dwarfs in the near-infrared and IRAC 
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Figure 10. The redshift probabihty distribution for Object 9697, 
which shows the best-fit photometric redshift of z = 6.92. 



bands. In order to rule out our candidates as brown dwarfs 
we compared their near-infrared colours to the following 
three colour spaces: [3.6 — 4.5]; [J — 4.5]; and [Ks — 4.5]. 
We combined the errors on our data points to conserva- 
tively explore the colour space covered by our objects and 
compared them with the expected brown dwarf colours. We 
found [3.6 — 4.5] colours indicative of spectral types T2-5 for 
object 9136 but we found [J — 4.5] colours indicative of an 
L8 dwarf and [Ks - 4.5] colours of a T7 or T8 dwarf. Each 
constraint contradicts the next, therefore we find it unlikely 
that object 9136 is a brown dwarf. 

Object 9697 has [3.6 — 4.5] colours consistent with all M, 
L and TO-3 spectral types and [Ks — 4.5] colours consistent 
with L5-8 and all T-dwarfs and (at the limit of its errors) 
[J — 4.5] colours between 3.47 and 4.59. The plots provided 
be Patten et al. (2006) indicate that a brown dwarf of type 
M, L or T will have [J — 4.5] colours of < 3.5, so at the 
extreme of its errors, object 9697 has colours just consistent 
with an L8 brown dwarf. However when we compared typical 
Y — J colours of brown dwarfs (Hewett et al. 2006) with 
Kg — 3.6/im colours, see Figure pT] we find that object 9136 
lies significantly away from the typical low-mass star colours, 
again making it unlikely that our objects can be explained 
as brown dwarfs. A combination of resolution (in the near- 
infrared bands) and signal-to- noise ratio (in the z^— band) is 
insufficient to determine whether this object is unresolved, 
as would be expected for a brown dwarf. 




Yab Jab 



Figure 11. y — Jvs. X — 3.6 colours for a sample of brown 
dwarfs from L1-T8. The points on the plot mark our three of 
our candidates, from left to right 9697, 2200 and 9136. Object 
9266 is not shown as it is only detected in the Y band. All 3 of 
these candidates lie away from the brown dwarfs marked by plus 
signs. Brown dwarf colours are taken from Leggett et al. (2000, 
2001,2002), Kirkpatrick et al. (2000), Geballe et al. (2002), Knapp 
et al. (2004) and Hewett et al. (2006). 



4.2 Plausibility of y-drop candidates 

Our final list of y-drop candidates consists of 4 objects. 
They span a magnitude range of Jab = 25.0 — 25.4 after 
applying aperture corrections. 

Object 2058 is detected only in the J-band with no IRAC 
source associated. 

Object 4551 is also only detected in the J-band but falls 
in a noisy region of the y-band image. 
Object 5512 falls on the edge of the ACS images. 
Object 4532 has a J-band detection which is possibly as- 
sociated with the extended edge of a galaxy 1^.^2 away. This 
source is detected in all of the ACS bands and has a match to 
GOODS-MUSIC object 9610 with a spectroscopic redshift of 
z = 0.347. It is likely that the large apparent (Y — J) colour 
recorded at the position of object 4532 is due to worse see- 
ing in the J-band than the y-band, or perhaps an intrinsic 
colour gradient in GOODS-MUSIC 9610. Given the proxim- 
ity of this low-redshift source, it is extremely unlikely that 
object 4532 is a genuine y-drop at 2; > 8. 

We find no robust y-drop candidates as all of our can- 
didates only appear in the J-band image with no significant 
Ks-hand detection which is unexpected as the J and Ks- 
bands probe similar depths. They also have no clear IRAC 
detection. This could indicate spurious detections in the J- 
band because they are unconfirmed in any other, or it could 
be the result of line contamination in the J-band filter from 
high equivalent width Lyman a. They could also be galaxies 
with very blue spectra indicating low metallicity and very 
little or no dust. 

In order to assess whether such blue colours could be 
plausible, we consider two scenarios, one of which has the 
colour difference produced by a blue spectral slope, the sec- 
ond assumes that the brightness in the J-band relative to Ks 
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is attributable to a very strong emission line (e.g. Lyman-a 
at 8.0 < z < 10.5). We then compare the constraints on 
the spectral slope and Lyman-a equivalent width with the 
known properties of Lyman-break galaxies at high redshift. 
If our limits fall outside the range observed in distant galax- 
ies, then the J-band detections are probably spurious (i.e., 
inconsistent with these being F-band drop-outs at 2; > 8). 

We assume a simple power law for the spectrum of 
fx oc (or equivalently fu oc A^^^, where /3 = — 2 is a 
spectrum of constant AB magnitude, flat in fi^). We place 
2cr upper limits on the spectral slope oi f3 < —2.44 (for 
sources 4551 and 5512 with Jab = 25.0) and P < -2.06 (for 
sources 2058 and 4532 with Jab = 25.4); these limits are 
conservative, because if the Lyman-a break occurs in the J- 
band filter (i.e. 2; > 8) rather than shortward of 1.1 fim then 
the true spectral slope would be even bluer. Our 2 a limits 
are consistent within the errors with the reported average 
for z = 6 i^-drop galaxies. Stanway, McMahon & Bunker 
(2005) derive values of (3 — — 2.2ili0.2 from i^-drop galaxies 
at z 6, with Bouwens et al. (2008) reporting (3 = —2.0. 

As an alternative to the blue {J — Ks) colour being due 
to a steep blue spectral slope, we now consider whether the 
apparent flux excess in the J-band could be due to emis- 
sion line contamination. To determine lower limits on the 
equivalent width of this putative line emission, we assume 
a spectrum flat in fi^ (i.e., constant AB magnitude with 
wavelength) longward of Lyman-a; this is typical of a low- 
extinction star-forming galaxy. We take our 2 a upper limit 
for the flux density {fu) in the K^-band as the upper limit 
on the continuum level in the J-band filter, and attribute 
the > 1.4x greater flux density (for objects 4551 and 5512) 
in the J-band (2cr lower limit) to be due to line emission. 
The J band has a width of 3000 A, which sets a 2 a lower 
limit on the observed equivalent width of EWobs > 1205 A 
(for J-drops 4551 & 5512). This corresponds to a rest-frame 
equivalent width of EWrest > 133 A if the line is Lyman-a 
at 2; = 8, where it enters the J-band (this limit is conser- 
vative because at larger redshifts, an increasing fraction of 
the J-band falls below the Lyman-a break, so the contri- 
bution of the continuum to the flux density would be even 
lower and hence the equivalent width higher). These equiv- 
alent widths are plausible for the Lyman-break and Lyman- 
a emitter populations at high redshift (e.g., Dawson et al. 
2004). We can rule out lower redshift lines with higher con- 
fldence (e.g.. Ho; at z ^ 0.9 would have an implausibly high 
EWrest > 634 A, and [O 11] 3727 A dit z ^ 2.3 would have 

EWrest > 365 A). 

We regard the marginal J-band detections to be highly 
suspect, although we cannot eliminate them from our selec- 
tion based on their blue spectral slopes or high equivalent 
width line emission falling in the J-band. They are uncon- 
firmed in any other band studied here and may be the re- 
sult of spurious detections. Deeper imaging in J, H and K 
and/or spectroscopy is required to confirm or disprove the 
nature of these candidates. 

4.3 Discussion of other work 

There are a few authors who have provided galaxy candi- 
dates at z > 6.5 using the Lyman Break technique. As men- 
tioned before, one of our 2;'-drop candidates, object 2200 
has also been identified by Stanway et al. (2008) but with 



a slight offset of 0^.^8. Two other objects were also identi- 
fied as possible high-redshift candidates in that paper, one 
of which we would not expect to detect because it does not 
lie within our field and the other object has a signal-to- 
noise ratio of less than our cut of S/N > 5 in the F-band. 
Within the limits of their study, Stanway et al. found the 
luminosity function at 2; = 7 to be consistent with predic- 
tions of the luminosity function from Bouwens et al. (2007) 
and McLure et al. (2009) at 2; ~ 6. Two other objects 
within the GOODS-South field were originally flagged by 
Mannucci et al. (2007) as high-redshift candidates but were 
later dismissed as brown dwarfs based on their morpholo- 
gies, Spitzer colours and spectroscopic information. These 
objects were included in our original catalogues as they do 
have the colours of high-redshift galaxies, and they were 
also identified in the GOODS-MUSIC catalogue with pho- 
tometric redshifts of z ^ 6.9 (objects 11002 & 7004 in the 
GOODS-MUSIC catalogue). Based on the non-detection of 
any credible candidates Mannucci et al. placed constraints 
on the UV luminosity function at z = 7 and claimed strong 
evolution in the luminosity function from 2: = 6 to 2: = 7. 
Further evidence for the evolution of the luminosity func- 
tion from 2; = 3.8to2; = 6^7is presented in Bouwens 
et al. (2004, 2005, 2008). No robust J-drops were presented 
but a number of 2;^-drop candidates were found. However 
the observations used in the Bouwens et al. (2008) study 
were much deeper than in other searches and some results 
even implied evolution from 2; = 6 to 2; = 7 and a potential 
luminosity function at 2; = 7 was derived while constraints 
were set on the luminosity function at 2; = 9 (the J-drop 
population). 

4.4 Implications for the UV Luminosity Function 

The number of robust candidates we detect can constrain 
the UV luminosity functions at z > 7. We compute the 
number of galaxies expected to be selected within our sur- 
vey area for several luminosity functions, derived from lower- 
redshift samples. A significant discrepancy between our ob- 
served number counts and those predicted would argue for 
strong evolution in the star-forming population with red- 
shift. 

To model the predicted number counts, we first adopt a 
simple model spectrum of a star-forming galaxy, where the 
rest-UV spectrum is approximately flat in longward of 
Lyman-a (i.e., f3 — —2 where /a oc A^, appropriate for star- 
forming galaxies at z ^ 6 - Stanway, McMahon & Bunker 
2005) and is severely attenuated below Lyman-o; due to the 
opacity of the intervening neutral hydrogen absorbers (we 
adopt an absorption of Da — 0.99 for z > 6.5). The Y- 
band filter is sensitive to the UV continuum longward of 
Lyman-a at 6.6 < z < 7.7, although at the higher redshifts 
the galaxies would have to be extremely luminous to appear 
in our magnitude-limited sample - not only is the luminos- 
ity distance greater, but also a smaller fraction of the filter 
bandpass lies above Lyman-a. We model this by considering 
small increments of redshift {Az = 0.1) between 2; = 6 and 
2: = 8, and for each redshift bin we calculate the number of 
galaxies expected as a function of limiting apparent magni- 
tude. We also determine the expected z —Y colours to assess 
whether our colour cut will select star-forming galaxies in 
that redshift bin. By summing over all the redshift bins, we 
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Figure 12. z^-band drop-out candidates. Each postage stamp is 5x5 arcsec in size. All candidates are undetected in the optical wavebands 
at the 2(7 level but have > 5cr detections in the F-band. Objects 9136 and 9697 are detected in all of the longer wavelength filters, excluding 
the Ks-hand for Object 9697, making these our most plaussible candidates. 



obtain the expected surface density of 2;^-drop galaxies as a 
function of magnitude. We then consider our exposure maps, 
and use this to compute the various areas of sky observed 
to different limiting magnitudes. For each area observed, we 
correct the predicted number counts for our measured com- 
pleteness as a function of magnitude (see Section [3.2[ ). By 
summing the expected number of galaxies above our colour 
cut {z' — Y > 1.0) and our significance threshold {S/N > 5 
in y-band) for each area of the survey, we obtained the to- 
tal number of 2;^-drop star-forming galaxies we would expect 
to find if the assumed luminosity function is appropriate at 
z^7. 



We can compare the number of expected galaxies de- 
rived from models with the number we actually detect. In 
this paper we will compare two UV luminosity functions, one 
derived by Steidel et al (1999) for Lyman-break galaxies at 
z = 3 (the [/-band dropouts), and the other by Bouwens 
et al. (2007) for the Lyman-break population at 2; ~ 6 
(the i^-band dropouts). The Steidel et al. UV luminosity 
function at z = 3 has a faint end slope of a = —1.6 and 
Lspji = 15.0 Moyr"^ and $ = 0.00138 Mpc"^, where L* is 
derived from the rest-frame UV around 1500 A and has been 
converted to an effective star formation rate using the rela- 
tion Luv = 8 X 10^'^ X SFR ergss"^ Hz~^ (Madau, Pozzetti 
& Dickinson 1998), appropriate for a Salpeter (1955) stellar 
initial mass function. At 2; ~ 6 the Bouwens et al. luminos- 
ity function shows strong evolution in L* from z ^ 3, with 
^SFR — 8.6M0yr~-^ (equivalent to 0.575 L^y at 2; = 3). 



The faint end slope is also steeper at 2; ^ 6 (a = —1.74) and 
$ = 0.001135 Mpc"^ (which is 0.82 0*^3). 

For the deepest region of our survey (pointing 1), we 
expect a 2;^-drop surface density brighter than our 50 per 
cent completeness limit (Yab < 25.9) of 0.373 arcmin"^ and 
0.066 arcmin"^ for the Steidel et al. (1999) and Bouwens 
et al. (2007) luminosity functions respectively. Accounting 
for completeness and the different depths as a function of 
survey area, the total numbers expected are 29.5 it 5.4 or 
5.2 dt 2.3 if the z = 3 or z = 6 luminosity functions are 
appropriate for the Lyman-break population at z ^ 7 (Fig- 
Clearly, as we have only two robust candidate z- 



14) 



drops (and at most three), we can strongly rule out a model 
where there is no evolution in the rest-frame UV luminosity 
function from z = 7 to z = 3, as the number of high-redshift 
galaxies is over-predicted by a factor of 10. There is some 
evidence for evolution from 2; = 6 to 2; = 7: it is likely that 
some or all of out z'-drop candidates are not at 2; = 7, and 
hence the observed number is at least a factor of two less 
than the prediction based on the z = 6 luminosity function, 
although the statistical significance of this is marginal given 
the small numbers. 

The same models were applied for the Y-drop candi- 
dates. Our 50 per cent completeness limit in the J-band is 
at Jab = 25.4, correcting for incompleteness means we ex- 
pect to find 10.5 ± 3.2 and 1.1 ± 1 F-drops for the Steidel 
et al. and Bouwens et al. luminosity functions respectively 
(Figure 15). All four of our candidates lie very close to our 



5cr J-band cut and object 4532 is possibly the edge of an 
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Figure 13. F-band drop out candidates. Each postage stamp is 5x5 arcsec in size. The candidates are undetected in the optical wavebands 
down to the 2cr level. Each F-drop candidate is detected in the J-band with S/N > 5 but have no clear detection in any of the longer 
wavelength bands. 



extended object ^ l'!2 away. They are also undetected in 
i^Tg-band which is of a similar depth to the J-band and hence 
the sources may be spurious or the result of line contamina- 
tion in the J-band. Although we cannot rule out all of our 
y-band candidates we do not believe them with a high de- 
gree of confidence. Thus we are again inconsistent with the 
Steidel et al. luminosity function implying evolution in the 
UV luminosity function between z = 8 and z — 3. Within 
the errors we are consistent with the Bouwens luminosity 
function at 2; = 6. 



At the limit of our survey for the z^-drop at z ~ 7 we 
begin to probe the z — Q Bouwens et al. luminosity function 
if there is little evolution, but in order to constrain the lu- 
minosity function more effectively, deeper and/or wider ob- 
servations are needed. This is a possibility with the VISTA 
Deep Extragalactic Observations (VIDEO) and the Ultra- 



VISTA surveys (Arnaboldi et al. 2007; see also Figures 14 
and 15). Due to its large area of ^ 12 sq. degrees, VIDEO 



will be able to probe the bright end of the luminosity func- 
tion and after five years will reach 5cr limiting magnitudes of 
Yab = 24.6 and Jab = 24.5. UltraVISTA goes considerably 
deeper (to Yab = 25.7 and Jab = 25.5 over 1.5 sq. degrees 
and Yab = 26.7 and Jab = 26.6 over 0.75 sq. degrees) but 
over a smaller area than VIDEO, and hence UltraVISTA 
will be more effective at measuring the position of the break 
and the slope of the faint end of the luminosity function (see 
Figure [M). 



5 CONCLUSION 

We have searched for high-redshift drop-out galaxies in the 
GOODS-South field using the new HAWK-I F-band data 
covering ^119 arcmin^. We have complemented this data 
with VLT ISAAC J and Ks images in addition to HST- 
ACS images in B, 1/, i and z' along with the deep Spitzer 
data in these fields. We employed a selection criteria of (Y — 
J)ab > 0.75 for the y-drops and {z' — Y)ab > 1-0 for the 
z^-drops and a S/N > 5. These catalogues were matched 
to the GOODS-MUSIC catalogue to eliminate objects with 
optical detections from our candidate lists. Each remaining 
candidate was inspected by eye to eliminate remaining data 
artifacts, spurious sources and optical detections. A total of 
4 y-drop candidates were found within our data. Due to 
the fact that all of our Y-diop candidates are very close 
to our S/N > 5 cut and are only significantly detected in 
the J-band, we do not believe any with a high degree of 
confidence. If none of our F-band sources are indeed real, 
then this demands significant evolution in the UV luminosity 
function since z = 3 based on the predictions by Steidel et 
al.(1999). 

We find 4 possible z' drop candidates, one of which. 
Object 2200, has a probable low redshift solution of Zphot — 
3.25 — 3.85. Another of our candidates. Object 9266, is only 
detected in the y-band. However we have 2 robust candi- 
dates. Objects 9136 and 9697, which have significant detec- 
tions in the IRAC wavebands and photometric redshifts of 
Zphot > 6.9. Our findings show evolution in the luminosity 
function since z — 3 but are, within the Poisson error, con- 
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Figure 14. Expected number of z^-drops, with the sohd hne as- 
suming the z = 3 luminosity function from Steidel et al. (1996) 
and the dashed line the z = 6 luminosity function from Bouwens 
et al. (2006). The dotted dashed line marks the phase space 
probed by the F-band data. The point denotes the two candi- 
dates we found. The dotted, hyphenated and solid lines mark 
the phase spaces that will be probed by the Ultra VISTA Deep, 
the Ultra VISTA shallow and the VISTA VIDEO surveys respec- 
tively. The luminosity functions shown here are not constrained 
at the bright end but with the depth and area of the VISTA sur- 
veys we will be able to measure the form of the function at these 
bright magnitudes. 
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Figure 15. Expected number of Y-drops, the various lines are 
the same as those presented in figure^] but using the J band as 
the long wavelength detection band. 



sistent with the Bouwens et al. (2007) z = 6 luminosity func- 
tion although small number statistics preclude from making 
any strong statements on the evolution in the Lyman-Break 
population at 2; > 6. In order to constrain the UV luminosity 
function at these high redshifts more effectively we need to 
search over a wider and/or deeper area, this will be possible 
with the combination of surveys such as VIDEO and Ultra 
VISTA (Arnaboldi et al. 2007). 
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